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Abstract 

Static and dynamical properties of weakly coupled antiferromagnetic spin 
chains are treated using a mean-field approximation for the interchain cou- 
pling and exact results for the resulting effective one-dimensional problem. 
Results for staggered magnetization, Neel temperature and spin wave exci- 
tations are in agreement with experiments on KCUF3. The existence of a 
narrow longitudinal mode is predicted. The results are in agreement with 

general scaling arguments, contrary to spin wave theory. 
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One-dimensional quantum spin chains are interesting objects to study for a number of 
reasons. On the one hand, experimental systems are generally very well described by simple 
yet nontrivial Hamiltonians involving very few unknown parameters. The standard example 
is the Heisenberg model, with only one free parameter, the exchange constant. Comparison 
between theory and experiment then becomes a particularly stringent test, as a variety of 
data have to be explained by one single parameter. Moreover, both a large number of exact 
theoretical results and powerful analytical and numerical methods are available, making this 
comparison particularly interesting. On the other hand, in spite of their simplicity, models of 
quantum spin chains have lead to a number of unexpected and unconventional predictions. 
For example, the exact solution of the antiferromagnetic spin-1/2 Heisenberg chain shows 
that the low-lying excitations are spin-1/2 objects [|TJ (now called spinons), quite different 
from standard spin waves. This prediction has been confirmed experimentally quite recently 
in KCUF3. Another example is Haldane's prediction of a gap in the excitation spectrum 
for integer-S 1 antiferromagnets H which again has found experimental confirmation. 

Strictly one-dimensional models of course do not exhibit phase transitions into states 
with a broken symmetry. It is nevertheless clear that any real compound, like KCuF 3 , [0,§|,[7|, 



Sr 2 Cu03, or Yb4As3 [ID] some form of interchain coupling is present. Then three- 
dimensional magnetic long-range order can appear below a Neel temperature T/y. In the 
present paper I show that in this case a conceptually simple approach, namely treating the 
interchain coupling in the mean-field approximation and treating the resulting effective one- 



dimensional problem as exactly as possible, |y],[T2| gives a coherent description of the ordered 
state and produces nontrivial quantitative predictions for static and dynamic quantities that 
can be successfully compared to experiments, in particular on KCUF3 where detailed neutron 
scattering results are available. 0.|J 

I start with the natural model for coupled spin-1/2 antiferromagnetic chains, namely a 
spatially anisotropic Heisenberg model for parallel chains forming a square lattice (i.e. the 
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lattice has tetragonal symmetry): 

H = J 2_j Si,T ' &i+l,r + J± $>i,r " -S^r+d • (1) 

i,v ^,r,(5 

Here i and r label lattice sites along the chain (z) and perpendicular (x, y) directions, 6 
is summed over the two nearest neighbor vectors in the transverse directions, and Si >r is a 
spin-1/2 operator at lattice site (i,v). The longitudinal and transverse exchange constants 
are J (> 0) and J± < 0, where in order to be close to the situation in KCUF3 a ferromagnetic 
interchain coupling is used (but all of the subsequent results apply with minor modifications 
also to J± > 0). 

I now treat antiferromagnetic order using a mean-field treatment of the interchain cou- 
pling. Assuming the order to be oriented along the ^-direction in spin space, the Hamiltonian 
(H) transforms into an effective single-chain problem described by 

H 1 = J^2S t - S i+ i -hJ2(-lYS* - 2NJ x m 2 . (2) 

i i 

Here iV is the number of sites in the chain, m = (— iy(S?) is the staggered magnetization, 
and h = —AJ±mQ. We thus have a one-dimensional antiferromagnet in an effective staggered 
field h, with the order parameter mo to be determined by minimizing the energy. The next 
step is to transform (||) into a fermionic model by a Jordan-Wigner transformation and to 
go to the continuum limit. The resulting fermionic model then is described by 

H 2 = J dz {-iv{^ L d x ^ L - ipid x iP R ) - h(^ R + ^Vl) + 2gipyy R iP L }-2NJ ± m 2 .(3) 

Here ipL,R are standard fermion field operators for left- and right-moving fermions. Eq.([3D 
can be recognized as the massive Thirring model for which an exact Bethe ansatz solution 
exists. [^3j The parameter v determines the velocity of the low-lying excitations, and by 
comparing with the exactly known (for h = 0) spectrum of Hi can be fixed as v — irJa/2 
where a is the lattice constant. The proper identification of g is more delicate because the 
underlying fermionic lattice model has rather strong interaction and a naive transition to the 
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continuum limit therefore is uncontrolled. One can however notice that the exact solution 



of H 2 has a mass gap [13 



A = he^ t / an7F \ , (4) 

7T( 7 -1)' K) 

where 71/7 = 2arccot(— g/2v) parameterizes the interaction, and A is the "rapidity cut- 
off". |0| Requiring further that the total particle number is independent of mo, one has 
A oc ln(l//i), i.e. A oc K 1 . On the other hand, standard scaling relations for the lattice 
Hamiltonian Hi imply a mass gap oc i.e. we have 7 = 2/3. Thus the parameters in H 2 
are 

g = 2v = tcJcl . (5) 

It is now straightforward to obtain the variation of the ground state energy per site with mo 
as 

E(m ) - E(m = 0) = -2J ± ml - ^r^ avn l /3 (^j > (6) 

where Uq is the fermion density, equal to 1/ (2a) in the original lattice model. Minimizing with 
respect to m one finds immediately the equilibrium value of the staggered magnetization 

as 

The mass gap then is 

A =^|j L |» 6.175| J ± \. (8) 

57T 

The results (|7p and (Q) will be compared to experiment below. 

We now turn to the spin dynamics. The appropriate generalization of the above mean- 
field approximation then is an RPA treatment of the interchain interaction. In the ordered 
phase translational symmetry is broken, and therefore there are umklapp processes coupling 
modes at wavevectors q and q + Q (where Q = (0, 0,7r/a)). The transverse susceptibility 
then is a 2 x 2 matrix given by 

X(9>W) = - —— ■ — r, 9 

1 - 2|J L |(cosg ;c + cos q y )x{Qz,^) 
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where 



(10) 



is the susceptibility matrix of the one-dimensional model H\. Explicit expressions for Xn.u 
are not known, however, a great deal can be learned from general symmetry properties of Hi 
and H 2 : (i) by spin rotational invariance the magnitude of the staggered magnetization is 
independent of the orientation of the staggered field in Hi; this implies Xnijt, 0) = — 1/(4 Jj_) 
and guarantees the existence of a Goldstone mode (spin wave) in the ordered state, (ii) From 
the equation of motions derived from H\ one finds Xn(0,c^) = {h 2 / uj 2 )[x n {^ , oj) — XnC^jO)], 
Xu(0, oS) = Xui^,^) = {h / oj)[xn{K , uj) — XnC^jO)], at q z = 0, 7r everything is thus determined 
by x„(7r,u;) alone, (iii) when all relevant energies (u;,A, ...) are much smaller than J the 
relativistic invariance of H 2 can be used: the operators S + (q z ~ it) and S + (q z ~ 0) are a 
Lorentz scalar and vector, respectively. Consequently, Xni^ + q, w) is a function of uj 2 — v 2 q 2 
only, and Xu(q,u) = {huj/{uj 2 — v 2 q 2 ))[x n ( 71 ' + q,&) — Xn( 7r , 0)]. The full three-dimensional 
x(q,ou) thus is entirely determined by Xn(^ + q,w) alone, as long as q z is in the vicinity of 
or 7r and u <^ J . 

We can further use the known spectrum of H 2 (13| to determine the form of Xni 71 ' + q, UJ )'- 
this function involves intermediate states where the ^-component of the magnetization has 
increased by AS* 2 = 1, corresponding to an added fermion in the Thirring model language. 
These excitation have energy A, leading to a pole. There are of course also multiple- 
excitation contributions, leading to a continuous spectrum. The lowest such excitation comes 
from a combination of the elementary AS Z = 1 excitation with an AS Z = excitation which 
also has energy A (see the discussion of the longitudinal excitations below), thus leading 
to a threshold at uj = 2 A. Notice that the elementary excitations creating this continuum 
are thus quite different from the spinons responsible for the lowest continuum of an isolated 
chain. p| Multiple-excitation continua of course also exist, with thresholds nA,n > 3. One 
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thus can write 

Xn(n + q,uj) = — — — - + f(uJ 2 -v 2 q 2 ) , (11) 

A z + v^q* — or 

where the unknown function f(x) contains contributions from the many-excitation continua, 
has a threshold singularity at x = 4A 2 , and is real below the threshold. The constant z 
ensures that Xn( 7r ? 0) = — 1/(4 J^). 

The excitation spectrum now is given by the singularities of x{Qi u )i ec L-(H)) an d in par- 
ticular the low-energy (a; < 2A) states are found from poles. One finds: (i) for propagation 
along the chain there is a spin wave mode with a>(0,0, q z ) = v(J±)\q z \, where the the spin 
wave velocity is only weakly affected by interchain coupling: e.g. neglecting / in eq. (|TTD (the 



"single mode approximation", SMA), v(J±) = v / Jl + h 2 / A 2 ; (ii) for transverse wavevector 
(it, 0) and more generally on the whole line cosq x + cosq y = the spin wave frequency is en- 
tirely determined by the one-dimensional result, and therefore from eq.flTTD uj(tc, 0, tt) = A, 
i.e. the mass gap is directly accessible experimentally and proportional to J±; (iii) in the 
SMA the spin-wave dispersion is uj(q) 2 = A 2 (l — (cosq x + cosq y )/2) + v 2 q 2 . In the relevant 
case | Jj_| J the qr-dependence of the transverse part that of standard spin wave theory. 
Taking into account / will lead to a modified relation in the transverse directions; (iv) in 
the SMA, the static susceptibility is x(0, 0) = 147r/(405J) ~ 0.109/ J. Taking into account 
/ will decrease this value. One can notice that for J± —>■ one expects to recover the purely 
one-dimensional result x(0, 0) = \/{ji 2 J) ~ 0.101/ J. This suggests that the SMA is a rather 
good description of low-energy properties. 

There also are longitudinal excitations, corresponding to oscillations of mo about its mean 
value. Within the RPA, the longitudinal susceptibility xl is given by a formula analogous 
to eq.(|9[). From eq.@ Xi( 7r ) 0) = — 1/(12J^), and the frequency dependence can again be 
obtained from the excitation spectrum of the massive Thirring model: here excitations with 
AS Z = intervene, i.e. particle-hole pairs. Naively, one expects a continuum above a gap 
2 A, however, because of the interaction in H2, the lowest excitation is actually an excitonic 
bound state at energy A, |L3[ i.e. at the same energy as the lowest AS Z = ±1 excitation. 



T4| The one-dimensional xl then has the same form as Xn, eg. (pTTj) , and in the SMA, there 
then is a longitudinal mode at u) 2 L (q) = A 2 [l — (cosq x + cosg y )/6] + v 2 q 2 . Going beyond the 
SMA, the already relatively weak dispersion in the transverse directions would be further 
reduced. Notice that at q± = (it, 0) longitudinal and transverse modes are degenerate. 

In this approximation, the longitudinal mode is well-defined. This is however an over- 
simplification: a longitudinal mode can decay into two spin waves. [|15|] To get a quantitative 
estimate of this effect I have obtained an effective Ginzburg-Landau description of the model 
via a Hubbard-Stratonovich transformation of the interchain interaction in eq. ([]]). At the 
Gaussian level expressions like eq . (|T0D for the transverse and longitudinal susceptibilities 
are obtained, but higher order corrections can also be studied systematically. In particular, 
a longitudinal excitation can now decay into two spinwaves, and there is also important 
spinwave-spinwave scattering. Spin rotation invariance imposes that the matrix elements 
for these processes at low energies are given by XL('K,0)/h and Xhij, 0)//i 2 , respectively. 
Taking these processes into account in an RPA-like fashion, the longitudinal susceptibility 
takes the form 

1 - 2|J J _|(cos<2 a . + cos q y ){xL{qz,uj) + E(g,wj) 
This form still retains one-dimensional relativistic invariance: Xl{(q±, Qz), — 
Xl((q±, 0), (u 2 — v 2 q 2 ) l l 2 ). Numerical results for Im[xi] (which determines the neutron 
scattering intensity) are shown in fig.|I|. One notices in particular a very sharp feature very 
close to the pole of the unrenormalized (S = 0) propagator. Numerically, the width of this 
peak is found of order 0.01A, the transverse mode remains thus rather well defined even 
when decay into spinwaves is taken into account. This can be attributed to the low frequency 
of the mode and the limited phase space available for decay into spinwaves. In addition to 
the peak, one also notices an incoherent background which for q z = n extends down to zero 
energy. 

Finally, the Neel temperature can be determined in this approach from the divergence of 
the static susceptibility at q = (0, 0, 7r) as a function of temperature. The single-chain sus- 



ceptibility is given by [|l6i|T7| x(7r, 0; T) = A/ (JT) ln 1//2 (AJ/T), where numerical calculations 



give A fh 0.32, A ps 5.8. Eq.@ then implies the relation 



|JJ = . (13) 

' ' 4Aln 1/2 (A J/T N ) V ; 

The above results for staggered magnetization, excitation spectrum, and Neel tempera- 
ture are now compared to experimental results on KCuF 3 . From the excitation spectrum 
above T/v the exchange constant along the chains is J = 34meV. From the measured |5| 
mo = 0.25 eq.(^) then gives J± = —J/16 = — 2.1meV. From eq.(|8|) then follows the first 
prediction of the present theory: A = 13meV. From eqs.(^) and eq. (|10|) A is the spin wave 
frequency at q x = (jr, 0,7r). Experimentally || uj(qi) = 11.5meV, in rather good agreement 
with the prediction. A second nontrivial prediction of the present theory is the existence of 



a longitudinal mode, at co>l(0, 0, it) y2/3A ps lOmeV. It is then tempting to associate the 
sharp rise in the energy-dependent neutron scattering intensity (which does not differenti- 
ate between transverse and longitudinal modes) observed around u = lOmeV |||7j with this 
mode. On the other hand, it is not clear how the two-spinwave process proposed in ref. 
which is rather featureless around uj = A (see fig.0), can account for this result. It would 
clearly be of interest to study this point in more detail. 

Using the estimated values for J, Jj_ eq.(|T3"D gives an estimate for the Neel temperature: 
T/v ~ 60K, overestimating the experimental result (TV = 39K) by about 50%. This discrep- 
ancy is in part due to the fact that logarithmic correction terms, enhancing the tendency 
to ordering, are included in eq. ([l~3|) , but neglected in the initial fermionic model, eq.(|3|). In 
the fermionic language, the logarithmic terms come from an extra umklapp interaction [19| 
duiipRipRipLipL + h.c.) in eq.@. These terms destroy the solvability of the fermionic model, 
but can be taken into account perturbatively. [I7j To lowest order then the h 4 ^ 3 term in 
eq.(||) is multiplied by a factor (1 + y ln(f /A)) 1 / 3 , with y = g u j{jxv). With y w 0.25 [T? 



the experimental value of mo leads to J± = — 1.6meV. From this A = 11.4meV, very close 
to the experimental value, and T/v = 47 K. The remaining discrepancy between theory and 
experimental values of T/v may well be related to thermal fluctuation effects neglected here 
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and which can be estimated to be on the level of 10%. [0 Experiment indeed indicates 
persistence of short-range order well above T N . [|,[7| 

It is also interesting to consider the case of Sr 2 Cu03. |8],|9| From the experimental Tjy ~ 
5K and the estimated J ~ 220meV one gets J± ~ 0.12meV. This then leads to a predicted 
spin wave energy at q = (7r,0,7r) of A = l.lmeV and m = 0.036. This last number is 
consistent with the experimental upper limit mo < 0.05. || It should however be noticed 
that in S^CuOs there is considerable structural anisotropy in the directions perpendicular 
to the chains, and the value derived for J± therefore is an average over transverse directions. 

The continuum limit used in the present treatment is valid if the correlation length 
£ = d/A ^ v/{jiT N ) is large compared to the lattice constant. For the two compounds 
discussed above one finds £ = 5a and £ = 320a, satisfying this criterion. The mean-field 



treatment of the interchain coupling leads to T N oc J±, m oc \/\J±\/J, up to logarithmic 



corrections. These relations are consistent with general scaling arguments. [20] On the 



other hand, the frequently used renormalized spin wave theory |2l| predicts T N oc J J\J±\, 
m oc l/ln(J/| J±\), in clear contradiction with scaling. Apart from the quantitatively 
satisfying description of KCuF 3 , the present theory thus also seems more satisfactory on 
grounds of consistency with generally valid arguments. 

I am grateful to A. Keren for interesting comments and to D. A. Tennant and A. Sandvik 
for stimulating correspondence. Laboratoire de Physique des Solides is a Laboratoire Associe 
au CNRS. 
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FIGURES 

FIG. 1. Frequency dependence of the imaginary part of the longitudinal susceptibility (in units 
of 1/| J]_|) at q = (0, 0, Tr+q z ) for q z = (full line), vq z = A/2 (dashed), vq z = A (dash-dotted), and 
vq z = 2 A (dash-double-dot). Also shown is the two-spinwave cross section at q z = 0, according 
to ref. Q (double-dash-dot). 
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